Solvable models for the gamma deformation having X(5)as limiting symmetry. 
Removing some drawbacks of the existent descriptions 
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Two solvable Hamiltonians for describing the dynamic gamma deformation, are proposed. The 
limiting case of each of them is the X(5) Hamiltonian. Analytical solutions for both energies and 
wave functions, which are periodic in 7, are presented in terms of spheroidal and Mathieu functions, 
respectively. Moreover, the gamma depending factors of the transition operator can be treated. 

PACS numbers: 21.60.Ev, 03.65. Ge, 21.10.Re 



Phenomenological formalisms like liquid drop (LD) , 
rotation- vibration Q , gamma unstable (GU) jj] , triaxial 
rotor (TR) HI Greiner-Gneuss (GG) [3, interacting bo- 
son (IBA) 0, 1 coherent state model y| , interpret some 
of the existent data by referring to a nuclear equilibrium 
shape, defining a nuclear phase with specific properties. 
It has been noticed that a given nuclear phase may be 
associated to a certain symmetry. Thus, the gamma un- 
stable nuclei can be described by the 0(6) symmetry 
the gamma triaxial nuclei by the rigid triaxial rotor D2 
symmetry [f|, the symmetric rotor by the SU(3) symme- 
try and the spherical vibrator by the U (5) symmetry. 

IBA succeeded to describe the basic properties of low 
lying states in a large number of nuclei in terms of the 
symmetries associated to the system of quadrupole (d) 
and monopolc (s) bosons which generate a U(6) algebra. 
The three limiting symmetries U(5), 0(6), SU(3) men- 
tioned above, are dynamic symmetries for U{6). More- 
over, for each of these symmetries a specific group reduc- 
tion chain provides the quantum numbers characterising 
the states, which are suitable for a certain region of nu- 
clei. A nice classification scheme was provided by Casten 
0, who placed all nuclei on the border of a symmetry 
triangle. The vertices of this triangle symbolise the U (5) 
(vibrator), 0(6) (gamma soft) and SU(3) (symmetric ro- 
tor), while the legs of the triangle denote the transitional 
region. In Ref . [TOl ITU , it has been proved that on the 
(7(5) — 0(6) transition leg there exists a critical point for 
a second order phase transition while the U (5) — SU (3) 
leg has a first order phase transition. 

Recently, Iachello El pointed out that these crit- 
ical points correspond to distinct symmetries, namely 
E(b) and X(b), respectively. Remarkable is the fact 
that the experimentalists found representatives for the 
two symmetries |l4l Il5j . Short after the pioneering pa- 
pers concerning critical point symmetries appeared, some 
other attempts have been performed, using other poten- 
tials like Coulomb, Kratzer [l^ and Davidson potentials 
[l7| . These potentials yield also Schrodinger solvable 
equations and the corresponding results may be inter- 
preted in terms of symmetry groups. 

The departure from the gamma unstable picture has 



been treated by several authors whose contributions are 
reviewed by Fortunato in Ref. 19]. The difficulty in 
treating the gamma degree of freedom consists in the fact 
that this variable is coupled to the rotation variables. A 
full solution for the Bohr-Mottelson Hamiltonian includ- 
ing an explicit treatment of gamma deformation variable 
can be found in Ref. [23, E| • Therein we treated sepa- 
rately also the gamma unstable and rotator Hamiltoni- 
ans. A more complete study of the rotor Hamiltonian and 
the distinct phases associate to a tilted moving rotator is 
given in Ref. (2^. In the recent publications, for the sake 
of simplicity, one uses model potentials which are sums 
of a beta and a gamma depending potentials. In this way 
the nice feature for the beta variable to be decoupled from 
the remaining 4 variables, specific to the harmonic liquid 
drop, is preserved. Further, the potential in gamma is 
expanded either around to 7 = or around 7 = 7r/6. In 
the first case, if only the singular term is retained one ob- 
tains the infinite square well model described by Bessel 
functions in gamma. If the 7 2 term is added to this term, 
the Laguerre functions are the eigenstates of the approx- 
imated gamma depending Hamiltonian, which results in 
defining the representation of the A (5) symmetry group. 
The drawback of these approximation consists in that 
the resulting function are not periodic as required by the 
starting Hamiltonian. Moreover, they are orthonormal- 
ized on unbound intervals although the underlying equa- 
tion was derived under the condition of I7I small. The 
scalar product for the space of the resulting functions is 
not defined with the measure | sin 37^7 as happens in the 
LD model. Under these circumstances the approximated 
Hamiltonian in 7 is not Hermitian. 

In the present letter we offer two solutions which cure 
the mentioned drawbacks. One is based on spheroidal 
functions while the other one on the Mathieu equation. 
These cases will be treated separately. However, in order 
to stress on the importance of curing the bad behaviour 
of the gamma depending wave function we start by say- 
ing few words about the the gamma expansion models. 
This helps us to fix the notations and introduce the basic 
equations needed here. Thus, let us consider the Hamil- 
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where [/ is a periodic function in 7 with the period equal 
to 2n and 
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with Qk denoting the components of the intrinsic an- 
gular momentum. Any approximation for the potential 
term, for example by expanding it in power series of 7, 
alters the periodic behaviour of the eigenfunction unless 
a special caution is taken. Moreover, the approximat- 
ing Hamiltonian loses its Hermiticity with respect to the 
scalar product defined with the measure for the gamma 
variable, | sin 37 1 dj. 

We illustrate this by considering the case of a little 
more complex potential 

u\ cos 37 -I- U2 cos 2 37. (3) 
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Substituting ip — y|sin 3y\ip, the eigenvalue equation 
Hip = Eip, becomes Hip = with 
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•C 1. Expanding the terms in 7, in 
power series up to the fourth order, one obtains: 
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The low index accompanying U and W suggests that 
the expansion was truncated at the fourth order. Note 
that due to the term W, the equations of motion for the 
variable 7 and Euler angles are coupled together. Such 
a coupling can in principle be handled as we did for the 
harmonic liquid drop in Ref. [2ClL l2ll| . Here, we separate 
the equation for 7 by averaging W4 with an eigenfunction 
for the intrinsic angular momentum squared. The final 
result for H4 is: 
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where L denotes the angular momentum. If the aver- 
age is made with the Wigner function D^ K , important 
simplifications are obtained since the following relations 
hold: 



(Ql-Ql) = 0, (Q 2 3 )=K 2 



(7) 



Let us stick to this situation. Note that H4 contains 
a singular term in 7 at, 7 = 0, coming from the term 
coupling the intrinsic variable 7 with the Euler angles. 
One could get rid of such a coupling term by starting 
with a potential in gamma containing a singular term 
which cancels the contribution produced by the W term. 
Thus, the new potential would be 
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The corresponding fourth order expansion of if is: 
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Some remarks concerning the equation H' A Lp = are 
worth to be mentioned. If in this equation one ignores 
the 7 4 term, the resulting equation has the Laguerre func- 
tions as solutions and moreover the Hamiltonian exhibits 
the X(5) symmetry. Note also that the Hamiltonian co- 
efficients are different from those of Refs.[13,19]. The 
difference is caused by the fact that here, the expan- 
sion W4 is complete. Taking in the expanded poten- 
tial u\ = ui — and ignoring, for 7 small, the term 
27if 2 7 2 /20, the resulting potential is that of an infinite 
square well which was treated by Iachello in Ref. [lflj . 
The solutions are, of course, the Bessel functions of half 
integer indices. None of the mentioned solutions is peri- 
odic. Also the approximated Hamiltonians are not Her- 
mitian in the Hilbert space of functions in gamma with 
the integration measure as introduced by LD. 

To overcome this principle problems, we try first to 
avoid making approximations. Thus, let us consider the 
Hamiltonian given by Eq. Q where instead of U we 
consider U' as defined by Eq. [S] and moreover ignore 
W. Changing the variable x — cos 37, the eigenvalue 
equation associated to this Hamiltonian becomes: 
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If u\ = U2 = K = 0, the solution is the Legendre poly- 
nomial P n while E = 9n(n + 1). Te case u\ = u 2 = 
was considered recently in Ref. |24j| . For other particular 
choice of the coefficients U\ and U2, the solution is readily 
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obtained if one compares the above equation with that 
characterising the spheroidal oblate functions |23 
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The prolate case is reached by changing c — > ic. For c = 
0, the solutions of Eq. 1)11(1 are the associated Legendre 
functions P™. For c ^ 0, S nm , with m,n integers and 
n > m > 0, are linear series of these functions. 

In the case u\ = 0, the solution of Eq. (10) is identi- 
fied as being the spheroidal function while the energy is 
related to X mn by a simple equation: 
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For |c| small the energies E nm exhibit the asymptotic 
expansion 
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Eq. I|13fl considered for a fixed m but various n, defines 
a band. Similar expansions may be derived for \c\ large. 
In order to save the space we give the 0(c -4 ) expansions 
for a few .E nTO 's. 
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Now, we shall focus on an approximate solution which 
preserves the periodicity in 7. To this aim, we consider 
the Hamiltonian 
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with A a positive constant. Changing the function by 
the transformation $ = I sin 

3 7 |-i/ 2 $ j for sin(37) ^ 0, 
the eigenvalue equation for H becomes: 
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Under the regime of I7I small, we take the 0(7 3 ) ex- 
pansion of the terms depending on 7 and in the final 
expression approximate 7 ~ sin 7. In this way, Eq. 1161 
becomes: 
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We suppose now that this equation is valid in the interval 
[0, 27r]. The equation (|17J) is just the trigonometric form 
of the spheroidal functions. The algebraic version is ob- 
tained by changing the variable x = cos 7. For A = 1, 
one obtains the Mathieu equation: 
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There are two sets of solutions, one even and one odd 
denoted by $ + (a, q,"f) and &~(a, q, 7), respectively. For 
5 = 0, both solutions are periodic, irrespective of a. 



$ + (a,0,7) = cos (Va7) , $ (a, 0, 7) = sin (Vai) 



(19) 
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FIG. 1: The characteristic curves c„ are plotted as functions 
of q for several values of n. 

For q ^ 0, the Mathieu functions are periodic in 7 only 
for a certain set of values of a, called characteristic val- 
ues. These are denoted by c+ for even and c~ for odd 
functions, respectively. In the plane (a, q), the character- 
istics curves c^ separate the stability regions, shown in 
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Fig.l by black colour, from the nonstability ones, indi- 
cated by white colour in the quoted figure. For q = the 



equalities c„(0) 



hold. By means of Eq. Q 7(1 the 



characteristic values determine the energy E. Thus, the 
energy spectrum is given by — u with = — ^-q. 
The corresponding wave functions are the elliptic cosine 
and elliptic sine functions respectively: 



Similar expansions can be derived also for \q\ 3> 1. 
Here we give only the energy expressions: 
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They form an orthogonal set. The matrix elements of the 
gamma depending factors of the transition operator can 
be easily calculated in Mathematica. Moreover, in the 
regime of |g|-small, these matrix elements can be analyt- 
ically calculated, since the following representation of the 
wave functions hold: 
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where n > 3,# + = and = 7r/2. The corresponding 
energies have the following expressions: 
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, n > 3. 



Obviously, a phase transition is determined by the 
combined effects coming from the behaviour of the wave 
function in the and 7 variables, respectively. For the 
X(5) symmetry, the (3 variable is described by a Bessel 
function of irrational index, while 7 by a Laguerre poly- 
nomial. Here, we propose to change the description of the 
7 variable either by a spheroidal or by a Mathieu func- 
tion. These functions are periodic and the corresponding 
Hamiltonians Hcrmitian. Moreover, in both versions, the 
X(5) Hamiltonian is obtained in the limit of small |7|. 
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